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ABSTRACT 



Aims. The chemistry of, and infrared (IR) emission from, polycyclic aromatic hydrocarbons (PAHs) in disks around Herbig Ae/Be 
and T Tauri stars are investigated. PAHs can exist in different charge states and they can bear different numbers of hydrogen atoms. 
The equilibrium (steady-state) distribution over all possible charge/hydrogenation states depends on the size and shape of the PAHs 
and on the physical properties of the star and surrounding disk. 

Methods. A chemistry model is created to calculate the equilibrium charge/hydrogenation distribution. Destruction of PAHs by 
ultraviolet (UV) photons, possibly in multi-photon absorption events, is taken into account. The chemistry model is coupled to a 
radiative transfer code to provide the physical parameters and to combine the PAH emission with the spectral energy distribution 
(SED) from the star+disk system. 

Results. Normally hydrogenated PAHs in Herbig Ae/Be disks account for most of the observed PAH emission, with neutral and 
positively ionized species contributing in roughly equal amounts. Close to the midplane, the PAHs are more strongly hydrogenated 
and negatively ionized, but these species do not contribute to the overall emission because of the low UV/optical flux deep inside 
the disk. PAHs of 50 carbon atoms are destroyed out to 100 AU in the disk's surface layer, and the resulting spatial extent of the 
emission does not agree well with observations. Rather, PAHs of about 100 carbon atoms or more are predicted to cause most of the 
observed emission. The emission is extended on a scale similar to that of the size of the disk, with the short-wavelength features less 
extended than the long-wavelength features. For similar wavelengths, the continuum emission is less extended than the PAH emission. 
Furthermore, the emission from T Tauri disks is much weaker and concentrated more towards the central star than that from Herbig 
Ae/Be disks. Positively ionized PAHs are predicted to be largely absent in T Tauri disks because of the weaker radiation field. 

Key words, astrochemistry - circumstellar matter - planetary systems: protoplanetary disks - infrared: general 



1. Introduction 



Polycyclic aromatic hyd rocarbons (PAHs; Leger & Puget 79841 
Allamandola et al. 1989) are ubiquitous in space and are seen in 
emission from a wide variety of sources, including the diffuse in- 
terstellar medium, photon-dominated r egions, circumstella r en- 
velopes, and (proto)planetary nebulae dPeeters et al.l 120041 and 
references therein). The PAHs in these sources are electronically 
excited by ultraviolet (UV) photons. Following internal conver- 
sion to a high vibrational level of the electronic ground state, 
they cool by emission in the C-H and C-C stretching and bend- 
ing modes at 3.3, 6.2, 7.7, 8.6, 11.3, 12.8 and 16.4 /im. 

Using the Infrared Space Observator y (ISO: iKessler et al] 

1996) , the Spitze r Space Telescope dWerner et al.l 12004a: 
Houck et al l [2004) and various ground-based telescopes, 
PAH features have also been obser ved in disks around 
Herbig Ae/Be and T Tauri stars J Van Kerckhove n et al 



2000; iHonv et all l200lf 
2003; van Bo ekel et al] 



iPeeters et all 120021: iPrzygodda et al 
2004bt lAcke & van den Anckerl 12004 



Geers et aT1l2006l) . Spatially resolved observations confirm that 



the emission comes from reg ions whose size is consistent with 
that of a circumstellar dis k (Ivan Boekel et al.ll2004al;lGeers et al.l 
120051: lHabart et alH 2006). Because of the optical and/or UV ra- 
diation required to excite the PAHs, their emission is thought to 
come mostly from the surface layers o f the disks dHabart et al.l 
2004b). A cke & van den Anckerl 12004) showed that PAH emis- 



sion is generally stro nger from flared disks (Me eus et al. 2001; 
iDominik et al.l l2003) than from flat or self-shadowed disks. 

Although the presence of such large molecules in disks and 
other astronomical environments is intrinsically interesting, it 
is also important to study PAHs for other reasons. They are a 
good diagnostic of the stellar radiation field and can be used 
to trace small dust particles in th e surface layers of d isks, both 
near the center and further out (Hab art et al.l 1200 4b). In addi- 
tion, they are strongly involved in the physical and chemical 
processes in disks. For instance, photoionization of PAHs pro- 
duces ener getic electrons, which are a major heating source 
of the gas dBakes & Tielensl ll994t iKamp & Dullemondl 120041: 
Jonk heid et al.ll2004l) . The absorption of UV radiation by PAHs 
in the surface layers influences radiation-driven processes closer 
to the midplane. Charge transfer of C + with neutral and neg- 
atively charged PA Hs affects the carbon chemistry. Finally, 
Hab art et al.l (l2004al) proposed PAHs as an important site of H2 
formation in photon-dominat ed regions, a process which is also 
important in disk chemistry (Jonk heid et al. 2006a b). Although 
this process is probably more efficient on grains and very large 
PAHs than it is on PAHs of up to 100 carbon atoms, the lat- 
ter may play an important role if the grains have grown to large 
sizes. 

Many PAHs and related species that are originally present 
in the parent molecular cloud, are able to survive the star 
formation process and eventually end up on planetary bod- 
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ies (TAllamand ola & Hudginsll2003l) . They add to the richness 
of the organochemi cal "broth" on planets in habitable zones 
dKasting et al.l Il993h . from which life may originate. Further 
enrichment is believed to come from the impact of comets. 
PAHs have now been detected in come tary material during the 
Deep Impact mission dLisse et al.l 120061) and r eturned to Earth 
by the Stardust mission (ISandford et al.ll2006f) . The icy grains 
that constitute comets also c ontain a variety of other molecules 
(Ehrenfreund & Fraser 2003). Radiation-induced chemical reac- 
tions between frozen-out PAHs and these molecules lead to a 
large variety of c omplex species, including some that are found 
in life on Earth (Berns tein et al.lll~999t lEhrenfreund & Se phton 
2006). These possibilities are another reason why it is important 
to study the presence and chemistry of PAHs in disks. 

of PAHs in an astronomical context 
with increasingly com plex an d ac- 
(lOmont 



The chemistry 
has been studied 
curate models in 



Lepp et al 
19961: ' 



1 r 



1988; 



the 

r 



last two decades 



Bakes & Tielensl 119941: ISalama et all 
Dartois & d'Hendecourj Il997b IVuong & Foingl \ 20Q(i 
IWeingartner & Draine! boOU 



1986; 



LePa geetal.1 1200 lL 120031: .... 

Bakes ctalj|2001albh : however, none of these were specifically 



targeted at PAHs in circumstellar disks. Disk chemistry models 
that do include PAHs only treat them in a ve ry simple manner 
(e.g. lJonkheid et al.ll2004l: Habart et al.ll2004bl) . In this paper, an 
extensive PAH chemistry model is coupled to a radiative t r ansfer 
model for circ umstellar disks (iDullemond & Dominikl [2004; 
iGeers et aT1l2006l Dullemond et al. in prep.). The chemistry part 
includes ionization (photoelectric emission), electron recombi- 
nation and attachment, photodissociation with loss of hydrogen 
and/or carbon, and hydrogen addition. Infrared (IR) emission 
from the PAHs is calculated taking multi-photon excitation into 
account, and added to the spectral energy distribution (SED) of 
the star+disk system. The model can in principle also be used 
to examine PAH chemistry and emission in other astronomical 
environments. 

We will present the chemistry model in Section [2] followed 
by a brief review of the radiative transfer model in Section [3] 
The results are discussed in Section [4] and our conclusions are 
summarized in Section[5] 



2. PAH model 

The chemistry part of ou r model is a comb ination of 
the models developed by iLe Page et all d2001l hereafter 
LPSB01) and IWeingartner & Drainei d2 00ll hereafter WD01 



Photodissociation is treated according to iLeger et al.l (Il989h . 
Where possible, theoretical rates are compared to recent ex- 
perimenta l data. Our model employs the new PAH cross sec- 
tions of Drai ne & Lil (120061, h ereafter DL06), which are an 
update of iLi & Draine d2001l) based on exper imental data 
(Matt ioda et al . 2005a b) and IR observations (e.g. lWerner et al.l 
2004b; Smit h et al. 2004). In this section, we present the main 
characteristics of our model. 



2.1. Characterization of PAHs 

The PAHs in our model are characterized by their number of 
carbon atoms, N c . A PAH bears the normal number of hydrogen 
atoms, N^, when one hydrogen atom is attached to each periph- 
eral carbon atom bonded to exactly two other carbon atoms (e.g. 



12 for coronene, C 24 H [2 ). The ratio between and A^ c is taken 
as (e.g. DL06): 



0.5 



f° = N° H /N C = \ 0-5 V257A^ 
0.25 



N c <25, 

25 <N C < 100, 

N c > 100 . 



(1) 



This formula produces values appropriate for compact (peri- 
condensed) PAHs. Elongated (catacondensed) PAHs will have a 
higher hydrogen coverage; however, they are believe d to be less 
stable and to convert into a more co mpact geometry (Wan g~et~ail 
ll997HDartois & d'Hendecourtll9"97l) . so only compact PAHs are 
assumed to be present. 

Every carbon atom that bears one hydrogen atom in the nor- 
mal case is assumed to be able to bear two in extreme conditions, 
so each PAH can exist in 2N^ + 1 possible hydrogenation states 
(0 < N n < 2N^). Furthermore, each PAH can exist in two or 
more charge states Z. The maximum and minimum attainable 
charge depend on the radiation field a nd the PAH's ionizatio n 
and autoionization potentials (WD01; Bakes & Tielensl fl994). 
The number of accessible charge states increases with PAH size. 

For a number of properties, the radius of the PAH is more im- 
portant than N c . The PAHs of interest are assumed to be spher- 
ically symmetric and, following e.g. WD01 and iDraine & Lil 
(2001), are assigned an effective radius a, 



0.468] 



1/3 



(2) 



where a is in A. This equation does not give the actual geometric 
radius; rather, it gives the radius of a pure graphite sphere con- 
taining the same number of carbon atoms. Equations (Q]) and (f2]) 
do not apply to very small PAHs (N c Sj 20), but like irregularly 
shaped PAHs, they are assumed not to be abundant enough to 
contribute to the emission (Section |2T6l >. 



2.2. Photoprocesses 

The absorption of a UV or visible photon of sufficient energy 
by a PAH causes either the emission of an electron (ionization 
or detachment) or a transition to an excited electronic and vi- 
brational state. In the second case, internal conversion and flu- 
orescence rapidly bring the molecule to a high vibrational level 
of the grou nd electronic state . From here, several processes can 
take place dLeger et alJll988l [19891) : (1) dissociation with loss 
of atomic or molecular hydrogen; (2) dissociation with loss of 
a carbon-bearing fragment; or (3) cooling by infrared emission. 
The cooli ng rate constant for a PAH with internal energy £jnt is 
given bv lLi&Luninel (120031) : 



J912A 



4;rfi, j(r[£ , int ])cr abs 
he I A 



dA, 



(3) 



with B,i(T[Ei nt ]) the Planck function at the PAH's vibrational 
temperature, T[E mt ] (IDraine & Lil 120011) . The cross sections, 
0"abs, will be treated in Section |231 

The yield of a single dissociation process i, T,, depends on 
the rate constants, k, of all possible processes: 



Yi = 



(4) 



To determine the rate F, of dissociation process 2, the yield is 
multiplied by the absorption rate, taking into account the possi- 
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bility of electron emission, and integrated over the relevant en- 
ergy range: 



r, 



I 

Jo 



(1 - Y em )Y iC r abs N ph dE 



(5) 



where F em is the photoelectric emission yield (Eq. (O) and N p h 
is the number of photons in units of cirT 2 s _1 erg -1 . The radi- 
ation field is treated explicitly at every point in the disk in an 
axisymmetric 3-D geometry, i.e., the wavelength dependence is 
taken into account as well as the magnitude. 

Photoelectric emission is the ejection of an electron from a 
PAH due to the absorption of a UV photon. The electron can 
come either from the valence band (photoionization; possible for 
all charge states) or from an energy level above the valence band 
(photodetachment; only for negatively charged PAHs). The pho- 
toelectric emission rate is given by WD01 as: 



= I TionO"absA^phd£' + I Y del o- del N ph dE , 



(6) 



where h is Planck's constant. Photoelectric emission will only 
occur when the photon energy exceeds the threshold of hv l0 n or 
hv det . The photodetachment yield, Fdet, is taken to be unity: every 
absorption of a photon with hv > hv det leads to the ejection of 
an electron. The photoionization yield, F; on , has a value between 
0.1 and 1 for hv ^ 8 eV and drops rapidly for lower energies. 
For photons with hv > hv lon , the photodetachment cross section, 
crjet, is about two orders of magnitude lower for N c = 50 than 
the ionization cross section, which is assumed equal to cr a t, s . 
In order to derive F e m as used in Eq. (0, Eq. © is rewritten: 



[ Y ion + T d e t (TzbsNpbdE . 

J0 \ 0"abs / 



(7) 



The quantity between brackets is the electron emission yield, 



cr det 

4 1 det > 

0"abs 



(8) 



taking Y\ w and Y det to be zero for photon energies less than hv\ oa 
and /ivdet, respectively. Furthermore, F em is not allowed to ex- 
ceed unity: each absorbed photon can eject only one electron. 

If no photoelectric emission takes place, the PAH can un- 
dergo dissociation with loss of carbon or hydrogen. Several 
theoretical schemes exist to calculate the dissociation rates. 
LPSB01 employed the Rice-Ramsperger-Kassel-Marcus quasi- 
equilibrium theory (RRKM-QET) and obtained hydrogen loss 
rates close to those determined experimentally for ben- 
zene (C fi H A nap hthalene (C 1() H g ) and anthracene (C 14 H 1Q ). 
Lege r et aiT(ll98 9) investigated the loss of carbon as well as hy- 
drogen, us ing an inverse Laplace transform of the Arrhenius law 
( Fors ffl972l) to determine the rates. In this method, the rate con- 
stant is zero when the PAH's internal energy, E' m y, originating 
from one or more UV photons, is less than the critical energy, 
Eq, for a particular loss channel. When E mt exceeds Eq, 



= A 



p(Ei„ t - £o,x) 

P(£int) 



(9) 



where Ax is the pre-exponential Arrhenius factor for channel 
X and p(E ) is the density of vibrational states at energy E. 
iLeger et alj dl989b considered the loss of H, C, C 2 and C 3 ; 
LPSB01 also took H 2 loss into a ccount for PAHs wi th N H > N^. 
Our model uses the method of Lege r et all dl989l) for all loss 
channels, with values for A and Eq given in TableQ] 



Table 1. Arrhenius parameters for the loss of carbon and hydro- 
gen fragments from PAHs. 3 



Fragment 




Z 


E (eV) 


A (s-') 


C 


all 


all 


7.37 


6.2 x 10 15 


c 2 


all 


all 


8.49 


3.5 x 10 17 


c 3 


all 


all 


7.97 


1.5 x 10 18 


H 


<N° H 


all 


4.65 


1.5 x 10 15 


H 


>N° H 


< 


1.1 


4 x 10 13 


H 


>N° W 


> 


2.8 


1 x 10 14 


H, 


>N° H 


<0 


1.5 


4x 10 13 


H 2 


>N° H 


>0 


3.1 


1 x 10 14 



E and A for C, C 2 and C 3 are from lLeger~ ;t al. ( 198 9]). E for H 
andH 2 is based on the RRKM-QET parameters from Le Pa ge et ail 
d200lh and modified slightly to obtain a better m atch with their 
rates. A for H is modified from Lege r et all dl989h and used also 
forH 2 . 



The intensity of the UV radiation field inside the disk is char- 
acterized in our model by G' Q = u vv /u^, where 



m uv = I iiydv = I 

J6 eV J(s eV 



13.6 eV 



(h/c)hvN p hdv 



(10) 



and = 5.33 x 10 14 erg cm 3 is the energ y density in the 



mean interstellar radiation field dHabing| |l968t). However, the 
radiation field inside the disk does not have the same spectral 
shape as the interstellar radiation field, so G' is used instead of 
G () to denote its integrated intensity. The exact UV field at ev- 
ery point and every wavelength is calculated by a Monte Carlo 
code, which follows the photons from the star into the disk in an 
axisymmetric 3D geometry (see also Section |3~TI ). As a result, 
regions of high optical depth may not receive enough photons 
between 6 and 13.6 eV to calculate G[ y In those cases, the lower 
limit of the integration range is extended in 1-eV steps and the 
energy density of the Habing field is recalculated accordingly. If 
G' {) still cannot be calculated between 1 and 13.6 eV, it is set to a 
value of 10~ 6 . 

In the strong radiation fields present in circumstellar disks 
(up to G' = 10 10 in the inner disk around a Herbig Ae/Be 
star), multiple photons are absorbed by a PAH before it can cool 
through emission of IR radiation. These multi-photon events re- 
sult in higher dissociation rates than given by Eqs. © and ([9j. 
For example, a PAH of 50 carbon atoms is photodestroyed in 
a radiation field of G' Q « 10 5 if multi-photon events are al- 
lowed, while G' ~ 10 14 is required in a pure single-photon 
treatment (Section lZBI ). However, the multi-photon treatment is 
much more computationally demanding, hence our model is lim- 
ited to single-photon treatments for all processes but photode- 
struction; the latter will be discussed in detail in Section |2"6l The 
errors introduced by not including a full multi-photon treatment 
will be discussed in Section l431 

As long as the internal energy of the PAH exceeds Eq for any 
of the loss channels, there is competition between dissociation 
and radiative stabilization. Typically, the emission of a single IR 
photon is not sufficient to bring E mt below Eq. If the emitted IR 
photons are assumed to have an average energy ^ir of 0.18 eV 
(6.9 //m; LPSB01), a total of n — {Ei at — #o)/?ir photons are 
required to stabilize the PAH. The /th p hoton in this cooling pro - 
cess is emitted at an approximate rate (Her bst & Dun bar 1991) 



^rad,/ — 73 



[ffjnt ~ Q" ~ 1)<?Ir] 

J0£ 



1.5 



(ID 
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with E mt in e V and s = 3(N C + N H ) - 6 the number of vibrational 
degrees of freedom. The competition between dissociation and 
IR emission occurs for every intermediate state. Henc e, the total 
radiative stabilization rate is dHerbst & Le PageF T999) 



&rad — k 



rad, 1 



n 

i=i 



k m d,i + &diss,H + ^diss.H, + &diss,C ' 



(12) 



where £di ss ,c denotes the sum of all three carbon loss channels. 
Note that Eqs. ( fTTT i and dl~2b are used only in the chemical part 
of the code, where the details of the cooling rate function are 
not important (LPSB01). The radiative transfer part employs the 
more accurate Eq. (|3). 

Alternative stabilization pathways such as inverse inter- 
nal conversion, inverse fluoresc ence and Poincare fluorescence 
dLeachlll987l: iLeger et al.H"l988l) are probably of minor impor- 
tance. Their effect can be approximated by increasing qm, thus 
creating an "effective" & ra( j that is somewhat larger than the "old" 
£ ra d. However, this will have no discernable effect on the rates 
from Eq. ©, so these alternative processes are ignored alto- 
gether. 

2.3. Absorption cross sections 

iLi & Draind (2001) performed a thorough examination of 
the absorption cross sections, <x abs , of PAHs across a large 
range of wavelength s. Following new exper imental data 
(Matt ioda et al.ll2005 a b) and IR observations ("e.g. lWerner et alJ 
2004b; Smith et al. 2004), an updated model was published in 
DL06, providing a set of equations that can readily be used in 
our model. The cross sections consist of a continuum contribu- 
tion that decreases towards longer wavelengths, superposed onto 
which are a number of Drude profiles to account for the cr-cr* 
transition at 72.2 nm, the 7r-7r* transition at 217.5 nm, the C-H 
stretching mode at 3.3 /mi, the C-C stretching modes at 6.2 and 
7.7 /mi, the C-H in-plane bending mode at 8.6 //m, and the C-H 
out-of-plane bending mode at 1 1 .2-1 1 .3 //m. Some of these pri- 
mary features are split into two or three subfeatures and several 
minor features are included in the 5-20 /mi range to give a bet- 
ter agreement with recent observations (DL06). The additional 
absorpt ion for ions in the near IR measured by Mattioda et al. 
(2005a,b) is included as a continuum term and three Drude pro- 
files. 

The absorption properties of a PAH depend on its charge. 
Neutral PAHs have 6.2, 7.7 and 8.6 /mi features that are a 
factor of a few weaker than do cations. The 3.3 /mi band 
strength increases from cati o ns to neut rals jBauschlicheri r2002; 
iHudgins et all 1200 lL l2000t ILanghoffl [1996J). The other fea- 
tures have similar intensities for each charge state. While the 
cation/neutral band ratios at 6.2, 7.7 and 8.6 /mi are mostly inde- 
pendent of size, those at 3.3 /mi decrease for larger PAHs. DL06 
provide integrated band strengths, cr. ,, for neutral and ionized 
PAHs, but they do not account for the size dependence of this 



feature in cations. Rather than taking cr 
DL06 do, our model uses 



z=\ 

int,3.3 



rr z=1 - rr z=0 x 

"inU^ ~ ^int,3.3 A 



1 + 



41 



- 14 



0-227^> 3 as 



(13) 



This way, the ionized/neutral band ratio approaches unity for 
large PAHs. The the oretical ratios for N c = 24, 54 and 96 from 
Bauschlicher (2002) are well reproduced by Eq. d!3t . 

LP SB01 based their opacities on experiments by iJoblinl 
d 19921) on PAHs in soot extracts. They fitted the cross section of 



an extract with an average PAH mass of 365 amu (correspond- 
ing to N c ~ 30) with a collection of Gaussian curves, obtaining 
practically the same values for A $ 0.25 /mi (E ^ 5 eV) as 
DL06. Between 0.25 and 0.5 /mi, the integrated cross section of 
LPSB01 is about 2.5 times stronger. Since LPSB01 were only 
interested in the UV opacities, their cross section throughout the 
rest of the visible and all of the infrared is zero. For our purposes, 
the different opacities in the 0.25-0.5 /mi range will only affect 
the hydrogen dissociation rates, and, as will be shown in Section 
14.51 the effects are negligible. 

The details of the infrared cross sections of negatively 
charged PAHs are not well under stood, with very different val- 
ues to be found in the literature (Bauschlicher & Bakes 2000; 
ILanghoffl [1996I) . However, this is not a problem in this work. 
Because of the small electron affinities, photodetachment is a 
very efficient process and absorption of a photon by a PAH anion 
leads to ejection of an electron rather than substantial emission in 
the infrared. Rapid electron attachment retrieves the original an- 
ion before the transient neutral PAH can absorb a second photon. 
Hence, if anions are present in steady state, they are excluded al- 
together when the infrared emission is calculated. 

Large PAHs can carry a double or triple positive charge in 
those regions of the disk where the optical/UV field is strong 
and the electron density is low. The differences between the cross 
sections of singly charged and neutral P AHs generally increase 
when going to multiply charged PAHs (Bauschlicher & Bakes 
2000; Ba kes et al.ll2001al) . However, no cross sections for multi- 
ply charged species exist that can be directly used in our model, 
so the cross sections for singly charged species are used instead. 
Since multiply charged species only constitute a very small part 
of the PAH population (Section 14. IK this approximation is not 
expected to cause major errors. 

2.4. Electron recombination and attachment 

Free electrons inside the disk can recombine with PAH cations 
or attach to PAH neutrals and anions; both processes will be re- 
ferred to as electron attachment. The electron attachment rate, 
r ea = a;ea«e, with n e the electron number density, depends on 
the frequency of collisions between electrons and PAHs and on 
the probability that a colliding electron sticks, as expressed by 
WD01: 



(8kT\ 

<*ea = Se 



1/2 



na J | t = r ,Z 



(14) 



where k is Boltzmann's constant, m e is the electron mass, so is 
the permittivity of vacuum, T is the gas temperature, e is the 
electron charge and t is the so-called reduced temperature. 

The sticking coefficient, s e , has a maximum value of 0.5 to 
allow for the possibility of elastic scattering. The electrons that 
do not scatter elastically, have to be retained by the PAH be- 
fore their momentum carries them back to infinity. The proba- 
bility of retention is approximately 1 - e~ a ^ 1 ", where Z e m 10 A 
can be considered the mean free path of the electron inside the 
molecule. WD01 include an additional factor 1/(1 + e 20 ~ N c) for 
PAH neutrals and anions to take into account the possibility that 
the PAH will be dissociated by the electron's excess energy, but 
this is only important for PAHs smaller than 24 carbon atoms 
(see also LPSB01) and they are assumed not to be present in the 
disk (Section l2T6l l. Thus, our model uses 

(15) 



s e = 0.5(1 -e- a/l ') 



for all PAHs with a charge Z > Z n 
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Expressions for J, a factor that takes into account the charge 
and size of the PAH and the temperature of the gas, can be found 
in lDraine & Sutinl (119871) . For neutral PAHs, J oz T ~ I/2 , so the 
attachment rate does not depend significantly on the temperature 
(Bak es &Tieieiisl [r994). J changes as T for PAH cations in 
the PAH size and gas temperature regimes of interest, leading to 
r ea oc r -1 / 2 . Finally, for PAH anions, both J and the attachment 
rate vary as e~ l t T . 

Recombination rates at 300 K for PAH cations of up to 
sixteen carbon atoms have been determined experimentally 
dBiennier et al.l2006tlNovotnv et al.l2005l:lHassouna et al.l2003t: 



Rebrion-Rowe et al.ll2003l: lAbouelaziz et al.ll 1993b . The experi- 
mental and theoretical rates agree to within a factor of ~ 2 for 
these small species, except for naphthalene (C 10 H„), where the 
theoretical value is larger by a factor of ~ 7. Experimental data 
on larger PAHs are needed to ascertain the accuracy of Eq. ([T4l > 
for the sizes used in our model. 

Experimental data on electron attachment to PAH neutrals 
are scarce. Rates for anthracene (C |4 H |f| ) have been rep orted 
from 9 x 10" 10 to 4.5 x 10~ 9 cm 3 s 1 dMoustefaoui et al.ll 19981 
Tobita et al.f l992). which do not depend greatly on temperature. 
Equation (TBI predicts a rate of ~ 2 x 10"'" cm 3 s _1 , about 
an order of magnitude lower, for temperatures between 10 and 
1000 K. For pyrene (C 16 H 10 ), however, the theoretical value, 
~ 2 x 10 -9 cm 3 s _1 , is an order of magnitude higher than the val- 
ues fo und experimentally, 2 and 4.2 x 10~ 10 cm 3 s _1 (Tobit a~et al.l 
11992b . Hence, we can only estimate that Eq. ( TPfl i is accurate to 
within about an order of magnitude. No experimental data are 
available on electron attachment to negatively charged PAHs or 
on electron recombination with PAHs carrying a multiple posi- 
tive charge, so the uncertainties in the theoretical rates for these 
reactions are at least as large. 

2.5. Hydrogen addition 

The addition of atomic hydrogen to neutral and ionized PAHs 
is an exothermic process requiring little or no activation en- 
ergy (Hirama et al. 2004; Bauschlicher 1998). The temperature- 
independent rates from LPSB01 for addition to cation s are used, 
which are based primarily on experimental work bv lSnow et alj 
(1998). Addition to neutrals is about two orders of magn itude 
slower for benzene dTriebertet al.ll 19981: iMebel et alii 1997b and 
has not been measured for larger PAHs. No rates are known for 
the addition of hydrogen to anions or to cations with Z > 1. We 

take *£tH = 1(r2 ^d = d,H and *2£h = k ldd,H- The rate depends on 
N H as described by LPSB01. 

It is assumed that molecular hydrogen can only attach to 
PAHs with N n < N^. Again, the rate for addition to cations 
from LPSB01 is used and divided by 100 for the neutral and 
anion rates. 

For naphthalene and larger PAHs, addition of atomic hydro- 
gen (PAH + H — > PAH + i) is assumed to be much faster than 
the bimolecular abstractio n channel (PAH + H — > PAH i + H2) 
dHerbst & Le Page] 1 1999b . so the latter is not included in our 
model. 



2.6. PAH growth and destruction 

The PAHs observed in disks are generally not believed to have 
been formed in situ after the collapse and main infall phases. 
Formation and growth of PAHs requires a hi gh temperature 
(~ 1000 K), density and acety l ene ab undance dCherchneff et alj 
Il992t iFrenklach & Feigelsonl 1 1989b and, especially for the 



smallest PAHs, a low UV radiation field to prevent rapid pho- 
todissociation. If these conditions exist at all in a disk, it is only 
in a thin slice (less than 0. 1 AU) right behind the inner rim, and 
it is unlikely that this affects the PAH population at larger radii. 
Another possib le method of late-stage PAH formation is in ac- 
cretion shocks dDesch & Connolly|2002b . but too little is known 
about these events to include them in the model. Hence, no in 
situ formation or growth of PAHs is assumed to take place. 

PAH destruction is governed by the loss of carbon fragments 
upon absorption of one or more UV photons. C-C bonds are a 
few eV stronger than C-H bonds (e.g. TableQ]), so carbon is lost 
only from completely dehydrogenated PAHs. In order to get an 
accurate destruction rate, multi-photon events have to be taken 
into a cco unt, as first recognized by iGuh athakurt a & Draind 
( 1989) and Sie benmorgen et al. dl992b . Our model uses a proce- 
dure based on lHabart et al.l d2004bh . A PAH ensemble residing 
in a given radiation field assumes a statistical distribution over 
a range of internal energies. This distribution is represented by 
P(£'i n t), which is normalized so that P(£ , i nt )d£ , j nt is the proba- 
bility to find the PAH in the energy interval Ei nt ..Ei Bt + dEi nt . At 
every E[ nt , there is competition between cooling and dissociation 
with loss of a carbon fragment. The probability for dissociation 
depends on the ratio between kc (the sum of the rate constants 
for all carbon loss channels in Eq. (0) and k\ R (the instantaneous 
IR emission rate, comparable to Eq. dTTT)): 



'7des(£'mt) 



kc(E ln t) 



(16) 



kc(Emt) + k lR (E int ) 

£ IR is a very flat function in the relevant energy regime, while kc 
is very steep, so 77d es is approximately a step function. 

The probability, /?des. that a PAH in a certain radiation field is 
destroyed is then found by integrating P(£ , int)?7des(£int) over a U 
energies: 



r 

Jo 



■P(£'int)'7de S (£ , int)d£ , i, 



(17) 



with P^int) calculated according to Guhathakurta & Draine 
d!989b . This is equal to the formula in lHabart et ail d2004bh if 
77de S is replaced by a step function and the internal energy is 
converted to a tem perature using the PAH's heat capacity (e.g. 
iDrain e & Li 2001). Destruction is assumed to take place if pd es 
exceeds a value of 10~ 8 . 

We define a critical radiation intensity, GJ, which is the in- 
tensity required to cause photodestruction of a given PAH within 
a typical disk lifetime of 3 Myr, i.e., the intensity required to get 
Tdiss.c = 1 /Tdiss.c = 3 Myr. The G* Q for both single-photon and 
multi-photon destruction is plotted in Fig. Q] Knowledge of the 
radiation field at every point of the disk (Fig. [3) is required to de- 
termine where PAHs of a certain size are destroyed. An approx- 
imate approach is to trace the radiation field along the r v ; s = 1 
surface, where the intensity will be shown to decrease almost as 
a power law (Section |4~2"1 ). In our Herbig Ae/Be model (Section 
13.21 ), PAHs of 50 carbon atoms (G* Q = 1.2 xlO 5 ) are destroyed out 
to 100 AU on the t v ; s = 1 surface. The destruction radius is larger 
for smaller PAHs and vice versa; e.g., PAHs with N c = 100 are 
only destroyed in the inner 5 AU. PAHs with less than 24 carbon 
atoms are not taken into account at all, because they are already 
destroyed when G' ~ 1 . Although there are regions inside the 
disk where the UV intensity is lower, the PAHs in such regions 
will not contribute significantly to the emission spectrum. For a 
T Tauri star, the radiation field in the disk is much weaker, so the 
destruction radius is smaller. In our model T Tauri disk, 50-C 
PAHs can survive everywhere but in the disk's inner 0.01 AU, 
while 100-C PAHs can survive even there (Section |4~6j >. 
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Fig. 1. Photodestruction of PAHs. (a) The solid line gives GJ, the 
radiation intensity for which the destruction timescale is shorter 
than the disk lifetime of 3 Myr, when multi-photon events are 
included. The dashed line gives G,*j when only single-photon 
destruction is allowed. The right vertical axis shows the corre- 
sponding "destruction radius" along the r v ; s = 1 surface in our 
model disk around a Herbig Ae/Be star (Section l4~2k PAHs are 
destroyed inwards of this radius. The disk's inner and outer ra- 
dius are indicated by the dotted lines. (b,c) The t v1s = 1 surface 
in a vertical cut through our model Herbig Ae/Be and T Tauri 
disks. The tick marks denote the radius inwards of which PAHs 
of a given size are destroyed within the disk lifetime; for the 
Herbig Ae/Be model, this corresponds to the right vertical axis 
of panel (a). PAHs of equal size can survive much closer to a T 
Tauri star than to a Herbig Ae/Be star. 



2. 7. Other chemical processes 

No reactions between PAHs and species other than H and H2 
are included in our model. Although the second-order rate coef- 
ficients for the addition of, e.g., atomic nitrogen and oxygen are 
comp arable to that for atomic hydrogen (LPSB01: ISnow et alJ 
119981) . the abundances of these heavier elements are not high 
enough to affe ct the chemical equ ilibrium. Formation of dimers 
and clusters (Rapa cioli et al.l 120051 an d references therein) and 
trapp ing of PAHs onto grains and ices (Gudipati & Allamandola 
120031) are also left out. The midplane of the disk, where den- 
sities are high enough and temperatures low enough for these 
processes to play a role, does not contribute significantly to the 
IR emission spectrum. 



3. Disk model 

3.1. Computational code 

The Monte Carlo radiati ve transfer code RADMC 
dDullemond & Dominikl |2004|) is used in combination with 
the more general code RADICAL dDullemond & Turollal feOQO) 
to produce the IR spectra from PAHs in circumstellar disks. 
Using an axisymmetric density structure, but following photons 
in all three dimensions, RADMC determines the dust temper- 
ature and radiation field at every point of the disk. RADICAL 
then calculates a spectrum from all or part of the disk, or an 
image at any given wavelength. 

The calculations in RADMC and RADICAL are based on 
the optical properties of a collection of carbon and silicate dust 
grains. Recently, PAHs were added as another type of grain to 
model the emission fro m the Herbig Ae star V V Ser and the 
surrounding nebulosity (Pontoppidan et al. 2006), f rom a sam- 
ple of Herbig Ae/Be and T Tauri stars ( Geers et"aT]|2006l) . and 
to study the effects of dust sedimentation (Dullemond et al. in 
prep.). The PAHs are excited in a quantized fashion b y UV and, 
to a lesser degree, visible photons dLi & Drainei r2002). and cool 
in a classical way according to the "conti nuous cooling" ap- 
pr oximation (|Guhafhakurta & Draindll989l) . which was found 
by iDraine & Lil d200ll) to be accurate even for small PAHs. A 
detailed description of t he PAH emission module is given in 
iPontoppidan et al. ( 20061). whe reas tests against other codes are 

described in iGeerset all 120061). 

In the models used by Pontopp idan et alJ d2006l) . lGeers et al.1 

(2006) and Dullemond et al. (in prep.), no PAH chemistry 
was included. PAHs of a given size existed in the same 
charge/hydrogenation state everywhere, or in a fixed ratio be- 
tween a limited number of states (e.g. 50% neutral, 50% ion- 
ized). In the current model, the chemistry is included in the fol- 
lowing way. First, a single charge/hydrogenation state for a given 
A^ c is included at a given abundance in the radiative transfer pro- 
cedure, to calculate the disk structure and radiation field at every 
point. Using these physical parameters, the equilibrium distribu- 
tion of the PAHs over all possible charge/hydrogenation states is 
then determined. After an optional second iteration of the radia- 
tive transfer to take into account the heating of thermal grains by 
emission from the PAHs, the spectrum or image is calculated. 

Some additional chemistry is added to the model in order 
to determine the electron and atomic and molecular hydrogen 
densities, which are needed to calculate the chemical equilib- 
rium of the PAHs. The electron abundance, x e = n e /rin, is set 
equal to the C + abundance, based on a simple equilibrium be- 
twe en the photoionization of neutral C and th e recombination of 
C + dBerginetalJ 120061 120031: iLeTeuffet ai]|2000l) . All hydro- 
gen is in atomic form at the edges of the disk and is converted 
to molecular form inside the disk as the amount of dissociating 
photons decreases due to self-shielding and shiel ding by dust 
(Ivan Zadelhoff et alj2003HDraine & Bertoldill996l). Th e H 2 for- 
mation rate is taken from lBlack & van Dishoeckl d 1987b . For the 
outer parts of the disk, which receive little radiation from the 
star, an interstellar radiation field with Go = 1 is included. 

The RADMC and RADICAL codes only treat isotropic scat- 
tering of£hotOTis ; _We_^ifodwith a different radiative transfer 
code (van Zadelhof f et al.ll2003l) that no significant changes oc- 
cur in the results when using a more realistic anisotropic scatter- 
ing function. 

Different PAHs (i.e., different sizes) can be included at the 
same time, and for each one, the equilibrium distribution over 
all charge/hydrogenation states is calculated. The abundance of 
each PAH with a given N c is equal throughout the disk, except in 
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those regions where G' {) > G^ (N r ) ; there, the abundance is set to 
zero. As discussed in Section |2~6l mixing processes are ignored. 

3.2. Template disk with PAHs 

For most of the calculations, a template Herbig Ae/Be star+disk 
model is used with the following parameters. The star has radius 
2.19Rq, mass 2.9\Mq and effective temperature 10 4 K, and its 
spectrum is described by a Kurucz model. No UV excess due to 
accretion or other processes is present. The mass of the disk is 
0.01 M©, with inner and outer radii of 0.48 and 300 AU. The in- 
ner radius corresponds to a dust evaporation temperature of 1700 
K. The disk is in vertical hydrostatic eq uilibrium, with a flaring 
shape and a slightly puffed-up inner rim (iDullemond & Dominikl 
2004). The dust temperature is calculated explicitly, whereas the 
gas temperature is put to a constant value of 300 K everywhere, 
appropriate for the upper layers from which most of the PAH 
emission originates. The results are not sensitive to the exact 
value of the gas temperature (see also Section l4~5l l. 

The same parameters are used for a template T Tauri 
star+disk model, except that the star has mass 0.58Mq and ef- 
fective temperature 4000 K. The dust evaporation temperature 
is kept at 1700 K, so the disk's inner radius moves inwards to 
0.077 AU. 

C 50 H 18 is present as a prototypical PAH, at a high abundance 
of 1 .6 x 10~ 6 PAH molecules per hydrogen nucleus to maximize 
the effects of changes in the model parameters. This abundance 
corresponds to 10% of the total dust mass and to 36% of the total 
amount of carbon in dust being locked up in this PAH, assuming 
an abunda nce of carbon in du st of 2.22 x 10~ 4 with respect to hy- 
drogen (Habart et al. 2004b). The model was also run for PAHs 
of 24 and 96 carbon atoms. 
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Fig. 2. Steady-state distribution of the most important 
charge/hydrogenation states of C 50 H V (N^ = 18) in a disk 
around a Herbig Ae/Be star (10 4 K). Each panel shows a cut 
through the disk, with the equator on the x axis and the pole on 
the z axis. The gray scale denotes the fraction at which a state is 
present compared to all possible states. The dashed line denotes 
the r v ; s = 1 surface. The dotted "emission line" connects the 
points where the PAH emission is strongest for a given distance 
R = Vx 2 + z 2 from the star. The thick black contour lines denote 
the region responsible for most of the PAH emission; from the 
outside inward, the contours contain 95, 90, 76, 52 and 28% 
of the total emitted power between 2.5 and 13.5 fim. They are 
omitted from the last frame for clarity. 



4. Results 

This work focuses on the chemistry of the PAHs in a circum- 
stellar disk and on its effects on the mid-IR emission, as well 
as on the diffe rences between d isks around Herbig Ae/Be and 
T Tauri stars. Geer s et ail (120061) analyzed the effects of chang- 
ing various disk parameters, such as the PAH abundance and the 
disk geometry; settling of PAHs and dust will be analyzed by 
Dullemond et al. (in prep.). 

4. 1 . PAH chemistry 

Due to the large variations in density and UV intensity through- 
out the disk, the PAHs are present in a large number of 
charge/hydrogenation states (Table |5J. When a disk containing 
only C 50 H V (A 7 ^ = 18) is in steady state, our model predicts that 
56% of the observed PAH emission between 2.5 and 13.5 fim 
originates from C 50 H lg , with another 9% from C 50 H|„. PAHs 
missing one hydrogen atom also contribute to the emission: 22% 
comes from C 50 H 17 and 12% from C 5Q H| 7 . 

A strong contribution to the observed emission does not im- 
ply a high abundance throughout the disk, because the contri- 
bution of each state to the spectrum also depends on its spatial 
distribution. Figure|2]shows a cut through the disk and indicates 
the abundance of the six most important states with respect to 
the total C 5() H y population. Near the surface, where the radiation 
field is strong, the PAHs are ionized and some of them have lost a 
hydrogen atom. Going to lower altitudes, the PAHs first become 
neutral and then negatively ionized. Still lower, the increasing 
density and optical depth lead to further hydrogenation, result- 



ing in a high abundance of the completely hydrogenated anion, 
C 50 Hr 6 , around the midplane. 

Since normally hydrogenated and positively ionized states 
occur in more optical/UV-intense regions than do completely hy- 
drogenated and negatively ionized states, the former will emit 
more strongly. For example, the state responsible for more than 
half of the emission, C^Hjg, constitutes only 21% of all PAHs in 
the entire disk (Table |2J. About half of the PAHs (45%) are ex- 
pected to be in the form of C sn H ^ fi . This state, which is assumed 
not to emit at all (Section I2.3K dominates the high-density re- 
gions close to the midplane. The normally hydrogenated anion, 
CjqHjj, accounts for 28% of all PAHs, and the remaining PAHs 
are mostly present as anions with 18 < N H < 36. States like 
C 5Q H n and C 50 H3 6 are entirely absent: the same environmental 
parameters that favour dehydrogenation (strong radiation, low 
density), also favour ionization. 

Species with N n < 17 are also predicted to be absent. They 
have to be formed from C 5Q H^ 7 , but wherever these states ex- 
ist, the ratio between the UV intensity and the hydrogen density 
favours hydrogen addition. In those regions where the ratio is 
favourable to photodissociation with H loss, the radiation field 
is also strong enough to destroy the carbon skeleton. The abun- 
dance of C 5Q H^ 7 is boosted by H2 loss from C 5Q H^ 9 , which is 
a much faster process than H loss from C 5() H^ g . Moving from 
the surface to the midplane, the increasing hh/Gq ratio allows 
all hydrogenation states from A^ + 1 to 2A^ to exist. However, 
since all of them except N n = 19 are negatively charged, they 
are assumed not to emit. 
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Fig. 3. The radiation field (characterized by G' Q ; Eq. ( [Tol l) for 
a flaring disk around a Herbig Ae/Be star (10 4 K), containing 
PAHs of 50 carbon atoms. The dashed line is the "emission line" 
from Fig. [2] The white contour lines are the same as the black 
contour lines in Fig. [2] 
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Fig. 4. Normalized cumulative PAH emission for the Herbig 
Ae/Be model as function of total hydrogen density (n H , left 
panel), radiation field (G' , middle panel) and electron density 
(n e , right panel). 



4.2. PAH emission 

The integrated UV intensity between 6 and 13.6 eV, character- 
ized by Gq (Eq. (fTUli), varies by many orders of magnitude from 
the disk's surface to the midplane. The radiation field for the 
model star+disk system is shown in Fig. [3] with the region re- 
sponsible for most of the 2.5-13.5 fim PAH emission indicated 
by contour lines. Part of this region is truncated sharply due 
to photodestruction of the PAHs, marking the G' = GJj line. 
An analysis of the disk structure shows that some 95% of the 
PAH emission originates from a region with 10 1 < G' Q < 10 5 , 
10 3 < n H (cm" 3 ) < 10 9 , and KT 1 < n e (cm 4 ) < 10 4 (see also 
Fig.gli. 

When the radiation field is traced along the T v j s = 1 surface, 
and the disk's inner 0.01 AU is disregarded, the integrated inten- 



sity decreases almost as a power law. Specifically, for any point 
on the T v is = 1 surface a distance R T (in AU) away from the star, 

G' {) ~ 3.0 x 10 8 /?; 1 ' 74 . (18) 

The power-law exponent differs from the value of -2 expected 
based on geometrical considerations because of the curvature in 
the T v ; s = 1 surface. This relationship allows one to predict to 
what radius R* T PAHs of a certain size are destroyed, as was done 
in Section [2~6l and Fig.fT] It should be noted, however, that this 
only applies to the T v j s = 1 surface. PAHs can survive and con- 
tribute to the emission from R < R* when they are at lower alti- 
tudes (e.g. Figs.[3]and|7]l. 

The abundances of the charge/hydrogenation states in the 
95% emission region are different than those in the entire disk. 
Although they agree more closely to the emissivity contribu- 
tions, a perfect correspondence is not achieved. For instance, 
only 1 .0% of the PAHs in this region are in the form of C 50 H| 8 , 
but they account for 8.0% of the emission (Tabled. This is again 
due to the positive ions occurring in a more intense radiation 
field than the neutrals and negative ions. 

In a vertical cut through the disk, a point exists for every 
distance R &m from the star where the PAH emission is strongest. 
These points form the "emission lines" in Figs.|2]and|3] The con- 
ditions along this line determine the charge and hydrogenation 
of the PAHs responsible for most of the observed emission. Its 
altitude depends mostly on two competing factors: the intensity 
of the UV field and the PAH density. A strong UV field leads to 
stronger emission, but a strong UV field can only exist in regions 
of low density, where the total emission is weaker. 

The variation of hh, n e and G' Q along the emission line is plot- 
ted in Fig. [5] In the disk's inner 100 AU (disregarding the actual 
inner rim), photodestruction of PAHs causes the emission line to 
lie below the t v - is = 1 surface, where «h and n e are relatively high 
(10 11 and 10 4 crrr 3 , respectively) and G' Q is relatively low (10 2 - 
10 3 ). The resulting n e /G' and hh/G,', ratios favour the neutral 
normally hydrogenated species, C 5Q H 18 . Going to larger radii, 
the emission line gradually moves up and crosses the t v ; s = 1 
surface at around 100 AU. As a consequence, «h and n e decrease 
while Gy increases. Because the electron abundance with respect 
to hydrogen also decreases, the PAHs first lose a hydrogen atom 
to become C 50 H 17 . Further out, ionization takes place to produce 
C 50 H| 7 . At still larger radii, the radiation field starts to lose in- 
tensity and C 5() H[ 8 becomes the dominant species. The deviation 
from the R~ 1 -' power law for the radiation field (Eq. ( fTBl ) is due 
to the emission line not following the r v j s = 1 surface. 

It should be noted that ionized and/or dehydrogenated 
species do exist (and emit) in the disk's inner part. However, 
for any given distance to the star, they only become the most 
abundant species at around 100 AU. 

4.3. Other PAHs 

If C 96 H 24 is put into the disk instead of C 5Q H 18 , the abundance 
of the dehydrogenated and normally hydrogenated states essen- 
tially goes to zero. This is due to the hydrogen addition rates be- 
ing larger and the hydrogen dissociation rates being smaller for 
larger PAHs. Most of the 2.5-13.5 fim emission in the N c = 96 
case comes from C % H 48 (44%), C 96 H| 8 (42%) and C 96 H 2 + 
(14%), while the anion, C % H 48 , is the most abundant overall 
(Table©. 

If only C 24 H 12 is put in, the PAH emission becomes very 
weak. The critical radiation intensity for this PAH is less than 
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Fig. 5. Total hydrogen density («hX electron density (n e ), radi- 
ation field (G' ) and ratio between electron density and radia- 
tion field (n e /G' Q ) along the emission line defined in Fig. |2]and 
Section l4~2l The bars at the bottom indicate the main emitter at 
each distance along the emission line. The parameters are those 
of the template Herbig Ae/Be model. 



unity (Fig. Q}, so it will only survive in strongly shielded ar- 
eas. There, 66% of the PAHs are present as C 24 H~ 2 , 28% as 
C 24 H 12 and 4% as C 24 H 24 . Only the neutral species contributes 
to the emission; however, since this PAH only emits from re- 
gions where G' {) ^ 1, no PAH features are visible in the cal- 
culated spectrum, despite the high abundance used in our model. 
This is exemplified in Fig. [6] where the calculated spectra for the 
model star+disk system are compared for the three PAH sizes. 
A disk around a Herbig Ae/Be star containing PAHs of 50 or 96 
carbon atoms shows strong PAH features, but the 24-C spectrum 
contains only thermal dust emission. 

The goal of Fig. [6] is to show the differences arising from 
the photochemical modelling of the three PAH sizes in model 
Herbig Ae/Be and T Tauri disks, rather than to provide realistic 
spectra from such objects. In order to fit our model results to 
observation s, one would need to include a range of PAH sizes in 
one m odel ( Li & Lunine 2003) a nd execute a larger parameter 
study dGeers et al.ll2006t Habart et al.ll2004bl) than was done in 
this work. 



4.4. Spatial extent of the PAH emission 

The top panel in Figure [7j presents the cumulative intensity of 
the five main PAH features and the continua at 3.1 and 19.6 
/mi as a function of radius for N c = 50. In accordance with 
Figs. |2] and [3] more than 95% of the power radiated in the fea- 
tures originates from outside the inner 10 AU, and some 80% 
from outside 100 AU. This is largely due to these PAHs being 
destroyed closer to the star. The 3.3 and 11.3 yum features are 
somewhat less extended than the other three features. The con- 
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Fig. 6. Model spectra (flux at 1 pc) for disks around a Herbig 
Ae/Be (black) and a T Tauri star (gray) containing either C 24 H [2 , 
C 50 H lg (shifted by a factor of 2) or C 96 H 24 (shifted by a factor 
of 4), with the charge and hydrogenation balance calculated by 
the chemistry model. 
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Fig. 7. Normalized cumulative integrated intensity of the five 
main PAH features (black) and the continua at 3.1 and 19.6 /mi 
(gray) for (top) C<- n H.„ around a Herbig Ae/Be star, (middle) 



'50^18 

C % H 24 around a Herbig Ae/Be star and (bottom) C 5Q H lg around 
a T Tauri star. In each case, the charge and hydrogenation bal- 
ance were calculated by the chemistry model. 



tinua are much more confined than the features, especially at 3.1 
fim, where 75% comes from the disk's inner rim. 

The PAH emission from a disk with N c = 96 (Fig. [7j mid- 
dle panel) is less extended than that from the same disk with 
N c = 50, because larger PAHs can survive at smaller radii. 
About 80% of the integrated intensity of the 3.3 /um feature orig- 
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inates from within 100 AU, and so does about 60% of the other, 
less energetic features. The continuum at 20 fim has the same 
spatial behaviour as the 3.3 yum PAH feature, while the contin- 
uum at 3.1 //m is very much confined towards the center. These 
results are in good ag reement with the spatial extent modelled by 
Hab art etaLI d2004bl) . The spatial extent of the continuum emis- 
sion from the C 96 H 24 disk is identical to that from the C 5() H Ig 
disk, so the thermal dust emission appears to be unaffected by 
the details of the PAHs and their chemistry. 

Our model also ag rees well with a number of spatially 
resolved observations dHabart et al.l 120061; iGeers et ail [2005; 
van Boekel et al. 2004a), where the PAH features were consis- 
tently found to be more extended than the adjacent continuum. 
Furthermore, PAH emission is typically observed on scales of 
tens of AUs, which cannot be well explained by our model if 
only PAHs of 50 carbon atoms are present. Hence, the observed 
emission is probably due to PAHs of at least about 100 carbon 
atoms. 

4.5. Sensitivity analysis 

The rates for most of the chemical reactions discussed in Section 
[2] are not yet well known. In order to gauge the importance of 
having an accurate rate for a given reaction, the equilibrium dis- 
tributions were calculated for rates increased or decreased by 
a factor of 100 from their normal model values. These results 
are also presented in Table [2] Modifying the hydrogen dissocia- 
tion and addition rates leads to no significant changes, so treat- 
ing photodissociation with loss of hydrogen in a purely single- 
photon fashion likely does not introduce large errors. The fact 
that most of the PAH emission comes from regions with a rel- 
atively weak UV field (G' Q < 10 5 ), where multi -photon events 
play only a minor role, further justifies the single-photon dehy- 
drogenation treatment. 

The ionization and electron attachment rates are more im- 
portant to know accurately. For instance, the contribution from 
the cations to the spectrum decreases by a factor of a few when 
taking an ionization rate that is 0.01 times the normal model 
rate. The shift away from positively charged species also affects 
hydrogenation (addition of hydrogen is faster to cations than 
to neutrals), resulting, e.g., in a smaller emissivity contribution 
from C 5Q H lg and C 5Q H[ g with respect to C 5Q H 17 and C 50 H| 7 . 
Further laboratory work on ionization and electron attachment 
rates, especially for larger PAHs, will help to better constrain 
this part of the model. 

As shown by e.g. iKamp & Dullemondl d2004l) and 
iJonkheid et al.l d2004l) . the gas temperature is not constant 
throughout the disk. The recombination rates between electrons 
and cations decrease for higher temperatures, whereas the 
attachment rates of electrons to neutrals are almost independent 
of temperature (Section |2~4l i. Taking a higher temperature would 
result in a slightly larger cation abundance. However, even for 
an extreme temperature of 5000 K, the recombination rates 
decrease by only a factor of ~ 4, and the effects will be smaller 
than what is depicted in Table|2] Hence, we believe it is justified 
to take a constant temperature of 300 K throughout the disk. 

4.6. TTauri stars 

The results discussed so far are for a star with an effective tem- 
perature of 10 4 K, appropriate for a Herbig Ae/Be type. Colder 
stars, like T Tauri types, are less efficient in inducing IR emis- 
sion in PAHs, as shown with both observations and models by 



IGeers et ail ((2006), unless they have excess UV over the stel- 
lar atmosphere. Hence, when the model parameters are other- 
wise unchanged, the observed flux from the star+disk system 
becomes lower and the absolute PAH features become weaker 
(Fig.©. 

It was shown in Section 14.21 that the integrated intensity of 
the UV field along the t v ; s = 1 surface of a disk around a Herbig 
Ae/Be star decreases approximately as a power law. This is also 
true for the T Tauri case, if the disk's inner 0.01 AU are again 
disregarded. The exponent is slightly larger: 

G' * 12R- Im , (19) 

with R in AU. The UV field peaks at G' Q = 6 x 10 5 at the inner 
rim and drops by two orders of magnitude within 0.01 AU, so a 
PAH of 50 carbon atoms (G* = 1.2 X 10 5 ) can survive almost 
everywhere (see also Fig.[TJi. 

The lower destruction radii also lead to the PAH emission 
being more concentrated towards the inner disk (Fig. [7] bottom 
panel). The 3.3 yum feature is particularly confined, with 70% 
originating from the inner 10 AU and 25% from the inner rim. 
As was the case for the Herbig Ae/Be disks (Section 14.4b . the 
inner rim contributes strongly to the continuum emission at 3.1 
jum. The 19.6 fim continuum also behaves very similarly to that 
from the Herbig Ae/Be disks, except that it is stretched inwards 
because of the smaller R in . Thus, the spatial extent of the thermal 
dust emission seems to be largely unaffected by the temperature 
of the central star. 

Cations are practically absent in the model T Tauri disk 
with any kind of PAH, due to the weaker radiation field. For 
N c = 24 and 50, almost all of the PAH emission originates from 
the normally hydrogenated neutral species, C 24 H [2 or C 5Q H lg . 
Less than 0.1% originates from other neutral states, primarily 
those missing one hydrogen atom (C 24 H n , C 5Q H 17 ) or having 
twice the normal number of hydrogen atoms (C 24 H 24 , C 5Q H 36 ). 
If a disk around a T Tauri star contains only 96-C PAHs, all 
of the emission is due to C 9g H 4g . The absence of cations could 
hel p explain the weak 7.7 and 8.6 /mi features in observed spec- 
tra dGeers et al.ll2006l) . because these two features are weaker in 
neutral PAHs. 

Anions are abundant for all three PAH sizes, accounting for 
about half of the entire PAH population. However, they are again 
assumed not to contribute to the emission. 

4.7. Comparison with observations 

lAcke & van den Anckerl d2004l) performed a comprehensive 
analysis of the PAH features in a large sample of Herbig Ae/Be 
stars observed with ISO, measuring line fluxes and compar- 
ing them to each other. Figure [8] recreates their Fig. 9, plot- 
ting the ratio of the integrated fluxes in the 8.6 and 6.2 fim 
bands against the ratio of the integrated fluxes in the 3.3 and 
6.2 fim bands. In order to gauge the plausibility of the numerous 
charge/hydrogenation states a PAH can in principle attain, the 
same ratios are also plotted for a sample of models containing 
one PAH in one specific state only. 

The 3.3/6.2 ratio from the model is very sensitive to the 
charge of the PAH and increases by an order of magnitude when 
going from ionized to neutral species. The 3.3 and 8.6 //m fea- 
tures are due to C-H vibrational modes, while the 6.2 fim fea- 
ture is due to a C-C mode (Section [ZBl , so both ratios in Fig. [8] 
increase with N H . The observations fall mostly in between the 
model points for neutral and ionized PAHs, so both charge states 
appear to contribute to the observed emission. This strength- 



Visser et al.: PAH chemistry and IR emission from circumstellar disks 



11 



2.0 




0.02 0.10 1.00 3.00 

P(3.3 ,am)/P(6.2 ^m) 



Fig. 8. The ratio of the fluxes of the 8.6 and 6.2 fim bands 
against the ratio of the fluxes of the 3.3 and 6.2 yum bands. 
Filled diamonds and arrows are detections and upper limits from 
lAcke & van den Anckerl d2004l) . Open diamonds and lines are 
predictions from the template Herbig Ae/Be model. The cross in 
the lower right shows the typical error bars for both the observa- 
tions and the model results. 

ens the model results presented in this work, although the ob- 
served emissivity contribution from neutral species seems to be 
somewhat less than the predicted 50%. Most of the observa- 
tions agree with the model prediction that the emission origi- 
nates from multiple hydrogenation states, with a lower limit of 
iV H = iV°-l. 

5. Conclusions 

The chemistry of, and infrared (IR) emission from, polycyclic 
aromatic hydrocarbons (PAHs) in disks around Herbig Ae/Be 
and T Tauri stars were studied. An extensive PAH chemistry 
model has been cr eate d, based primarily on the m odels of 
iLe Page et ail d200ll) and lWeingartner & Draiiiel d200ll) . with ab- 
sorption cross sections from lDraine & Lil d2006l) . This model in- 
cludes reactions affecting the charge (ionization, electron recom- 
bination, electron attachment) and hydrogen coverage (photodis- 
sociation with hydrogen loss, hydrogen addition) of PAHs in an 
astronomical environment. Destruction of PAHs by UV radiation 
is also taken into account, including destruction by multi -photon 
absorption events. By coupling the chemistry model to an exist- 
ing radiative transfer model, equilibrium charge/hydrogenation 
distributions throughout the disks were obtained. The main re- 
sults are as follows: 

• Very small PAHs (24 carbon atoms) are destroyed within a 
typical disk lifetime of 3 Myr even in regions of low UV 
intensity (G' ~ 1). No features are seen in the calculated 
spectrum for either a Herbig Ae/Be or a T Tauri disk, despite 
a high PAH abundance. 

• PAHs of intermediate size (50 carbon atoms) do produce 
clearly visible features, even though they are still photode- 
stroyed out to about 100 AU in the surface layers of a disk 
around a Herbig Ae/Be star. The model predicts that most of 
the emission arises from the surface layers and from large 
radii (more than 100 AU). Neutral and positively ionized 
species, bearing the normal number of hydrogen atoms or 
one less, contribute in roughly equal amounts. Negatively 



charged species are also present, but are assumed not to con- 
tribute to the emission. 

• Going to still larger PAHs (96 carbon atoms), photodestruc- 
tion becomes a slower process and the PAHs can survive 
down to 5 AU from a Herbig Ae/Be star. The slower pho- 
todissociation rates also mean that these PAHs are fully 
hydrogenated everywhere in the disk. Neutral and ionized 
species still contribute in comparable amounts to the emis- 
sion, with some 15% originating from doubly ionized PAHs. 

• The PAH emission is predicted to be extended on a scale sim- 
ilar to the size of the disk, with the features at longer wave- 
lengths contributing more in the outer parts and the features 
at shorter wavelengths contributing more in the inner parts. 
For similar wavelengths, the continuum emission is less ex- 
tended than the PAH emission. 

• Disks around T Tauri stars show weaker PAH features than 
do disks around Herbig Ae/Be stars because of the weaker 
radiation from T Tauri stars, assuming they have no excess 
UV over the stellar atmosphere. The PAH emission from T 
Tauri disks is considerably more confined towards the cen- 
ter than that from disks around Herbig Ae/Be stars, because 
PAHs can survive much closer to the star. For instance, a 
50-C PAH survives everywhere but in the disk's innermost 
0.01 AU. Furthermore, the radiation field is no longer strong 
enough to ionize the PAHs, and all the PAH emission origi- 
nates from neutral species for all three PAH sizes. This could 
help explain the weak 7.7 and 8.6 //m features in observed 
spectra. About half of all PAHs in a T Tauri disk are pre- 
dicted to be negatively ionized. 

• Comp aring the model results to spatially resolved observa- 
tions dHabart et alj2006llGeers et al.ll2005llvan Boekel et al.l 
2004a|) for Herbig Ae/Be stars, it appears that PAHs 
of at least about 100 carbon atoms are responsible 
for most of the emission. The emission from smaller 
species is predicted to be too e xtended. Other observations 
dAcke & van den Ancker! 2004) support the conclusion that 
the emission is due to a mix of neutral and singly positively 
ionized species. 
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Table 2. Abundances (columns 3 to 7) and 2.5-13.5 emissivity contributions (columns 8 to 12) of the dominant charge/hydrogenation states 
for three PAHs in the model Herbig Ae/Be star+disk system. ab 



Abundance in entire disk (%) c Contribution to 2.5-13.5 jim emission (%) c 





Z 


N.R. 


ioor ea 


o.oir ea 


ioor diss , H 


o.oir diss , H 


N.R. 


ioor ea 


o.oir ea 


ioor diss , H 


o.oir diss , H 






o.oir em 


ioor em 


o.oir add , H 


ioor add , H 




o.oir em 


ioor em 


o.oir add>H 


ioor add , H 












N r = 24, N° = 12, - 


1 < Z < 42 










12 


-1 


6.6(+l) 


9.2(+l) 


3.5(+0) 


6.7(4-1) 


6.5(41) 


- 


- 


- 


- 


- 


12 





2.8(+l) 


6.0(-l) 


8.9(+l) 


2.7(4-1) 
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1.0(+2) 


1.0(+2) 


1.0(42) 


1.0(42) 


1.0(42) 


13 


-1 


l.l(-l) 
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- 


14 


-1 


l.O(-l) 


2.4(-3) 
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3.0(40) 


6.8(40) 


1.2(41) 


17 


+2 


4.3(-4) 


8.7(-6) 


9.7(-4) 


1.5(-4) 


6.1(-4) 


1.3(-1) 


4.9(-3) 


2.4(-l) 


2.4(-2) 


1.9(-1) 


18 


-1 


2.8(+l) 


4.4(+l) 


1.1 (+0) 


3.0(41) 


2.6(41) 


- 


- 


- 


- 


- 


18 





2.1(+1) 


6.2(-l) 


4.8(+l) 


2.1(4-1) 


2.0(41) 


5.6(4-1) 


4.6(41) 


2.9(41) 


5.9(41) 


5.5(41) 


18 


+1 


3.3(-2) 


1.3(-3) 


3.7(-l) 


4.9(-2) 


3.1(-2) 


9.0(40) 


1.2(40) 


1.0(41) 


1.5(41) 


8.3(40) 


18 


+2 


7.5(-4) 


8.9(-6) 


8.3(-3) 


1.1 (-3) 


4.3(-4) 


2.3(-l) 


5.4(-3) 


2.4(40) 


3.51-1) 


1.3(-1) 


19 


-1 


9.6(-l) 


l.l(+0) 


7.6(-2) 


3.2(-l) 


5.6(-l) 


- 


- 


- 


- 


- 


19 





1.8(-1) 


3.3(-3) 


l.l(+0) 


5.3(-2) 


2.0(-l) 


5.3(-6) 


6.9(-6) 


9.K-6) 


3.9(-7) 


8.2(-l) 


19 


+ 1 


5.6(-3) 


3.0(-6) 


4.7(-l) 


3.1(-3) 


7.5(-3) 


1.7(40) 


- 


4.5(41) 


8.1(-1) 


2.4(40) 


20 


-1 


1.9(-1) 


1.4(-1) 


2.0(-3) 


1.2(-1) 


1.4(-1) 


- 


- 


- 


- 


- 


21 


-1 


1.8(-D 


5.2(-2) 


1.7(-3) 


1.2(-1) 


1.3(-1) 


- 


- 


- 


- 


- 


22 


-1 


i.sc-i) 


2.5(-2) 


1.7(-3) 


1.2(-1) 


1.3(-1) 


- 


- 


- 


- 


- 


23 


-1 


1.8(-1) 


1.6(-2) 


1.8(-3) 


1.2(-1) 


1.4(-1) 


- 


- 


- 


- 


- 


24 


-1 


1.9(-1) 


1.3(-2) 


1.8(-3) 


1.2(-1) 


1.4(-1) 


- 


- 


- 


- 


- 


25 


-1 


2.0(-l) 


1.2(-2) 


1.8(-3) 


1.2(-1) 


1.4(-1) 


- 


- 


- 


- 


- 


26 


-1 


2.0(-l) 


1.1 (-2) 


1.9(-3) 


1.2(-1) 


1.4(-1) 


- 


- 


- 


- 


- 


27 


-1 


2.1(-1) 


1.1 (-2) 


1.9(-3) 


1.3(-1) 


1.4(-1) 


- 


- 


- 


- 


- 


28 


-1 


2.2(-l) 


l.l(-2) 


2.0(-3) 


1.3(-1) 


1.5(-1) 


- 


- 


- 


- 


- 


29 


-1 


2.2(-l) 


1.2(-2) 


2.0(-3) 


1.3(-1) 


1.5(-1) 


- 


- 


- 


- 


- 


30 


-1 


2.3(-l) 


1.3(-2) 


2.1(-3) 


1.3(-1) 


1.5(-1) 


- 


- 


- 


- 


- 


31 


-1 


2.4(-l) 


1.4(-2) 


2.2(-3) 


1.3(-1) 


1.6(-1) 


- 


- 


- 


- 


- 


32 


-1 


2.5(-l) 


1.6(-2) 


2.2(-3) 


1.4(-1) 


1.7(-1) 


- 


- 


- 


- 


- 


33 


-1 


2.6(-l) 


2.7(-2) 


2.3(-3) 


1.4(-1) 


1.8(-1) 


- 


- 


- 


- 


- 


34 


-1 


2.7(-l) 


8.3(-2) 


2.4(-3) 


1.4(-1) 


2.2(-l) 


- 


- 


- 


- 


- 


35 


-1 


2.8(-l) 


4.6(-l) 


2.4(-3) 


1.4(~1) 


5.0(-l) 












36 


-1 


4.5(+l) 


5.3(+l) 


4.2(+l) 


4.5(4-1) 


4.8(41) 












36 





4.0(-l) 


3.0(-2) 


5.6(+0) 


3.8(-l) 


8.3(-l) 


5.0(-5) 


3.5(-5) 


6.2(-5) 


5.0(-5) 


5.2(-5) 












N r = 96, N° = 24, - 


2 < Z < 44 










48 


-2 


2.1(+1) 


2.0(+l) 


2.3(+l) 


2.1(4-1) 


2.1(41) 












48 


-1 


5.7(+l) 


7.9(+l) 


2.4(+l) 


5.7(4-1) 


5.7(41) 












48 





2.2(+l) 


7.1(-1) 


5.0(+l) 


2.2(4-1) 


2.2(41) 


4.4(41) 


8.7(41) 


1.3(41) 


4.4(41) 


4.4(41) 


48 


+ 1 


1.3(-1) 


1.0(-2) 


l.l(+0) 


1.3(-1) 


1.3(-1) 


4.2(41) 


1.2(41) 


2.9(41) 


4.2(41) 


4.2(41) 


48 


+2 


2.0(-2) 


3.8(-4) 


2.1(-1) 


2.0(-2) 


2.0(-2) 


1.4(41) 


5.0(-l) 


5.6(41) 


1.4(41) 


1.4(41) 



a N.R.: normal rates, i.e., the rates as discussed in Section|2] 100rx and O.OlFx: the normal rate for process X in/decreased by 
a factor of 100. The effect of increasing a certain rate is the same as decreasing the rate of the reverse process, ea: electron 
attachment; em: photoelectric emission; diss,H: photodissociation with loss of H or H2; add,H: addition of H or H2. 

b Only those states are shown whose abundance or emissivity contribution is at least 0.1% for the normal rates. 

c Percentage fraction of the entire PAH population in the disk. 



